Non-adiabatic Kohn-anomaly in a doped graphene monolayer 
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We compute, from first-principles, the frequency of the £2^, T phonon (Raman G band) of 
graphene, as a function of the charge doping. Calculations are done using i) the adiabatic Born- 
Oppenheimer approximation and ii) time-dependent perturbation theory to explore dynamic efi'ects 
beyond this approximation. The two approaches provide very different results. While, the adiabatic 
phonon frequency weakly depends on the doping, the dynamic one rapidly varies because of a Kohn 
anomaly. The adiabatic approximation is considered valid in most materials. Here, we show that 
doped graphene is a spectacular example where this approximation miserably fails. 
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Graphene is a 2-dimensional plane of carbon atoms 
arranged in a honeycomb lattice. The recent demonstra- 
tion of a field-effect transistor (FET) based on a few- 
layers graphene sheet has boosted the interest in this 
system [l, 2, 3]. In particular, by tuning the FET gate- 
voltage Vg it is possible to dope graphene by adding 
an excess surface electron charge. The actual possibil- 
ity of building a FET with just one graphene monolayer 
maximizes the excess charge corresponding to a single 
atom in the sheet. In a FET-based experiment, graphene 
can be doped up to 3 10^'^ cm~^ electron concentra- 
tion [l|, corresponding, in a monolayer, to a 0.2% 
valence charge variation. The resulting chemical-bond 
modification could induce a variation of bond- lengths and 
phonon-frequencies of the same order, which would be 
measurable. This would realize the dream of tuning the 
chemistry, within an electronic device, by varying Vg. 

The presence of Kohn anomalies (KAs) 0, Q in 
graphene could act as a magnifying glass, leading to a 
variation of the optical phonon-frequencies much larger 
than the 0.2% expected in conventional systems. On 
the other hand, the phonon-frequency change induced 
by FET-doping could provide a much more precise de- 
termination of the KA, with respect to other experimen- 
tal settings. KAs manifest as a sudden change in the 
phonon dispersion for a wavevector q ~ 2k^, where \ip 
is a Fermi-surface wavevector The KA can be deter- 
mined by studying the phonon frequency as a function of 
q by, e.g., inelastic x-ray, or neutron scattering. These 
techniques have a finite resolution, in q and energy, which 
limits the precision on the measured KA dispersion. In 
graphene, 2ki? is proportional to Vg. This suggests an 
alternative way to study the KA, that is to measure the 
phonon frequency at a fixed q and to vary 2k p by chang- 
ing Vg . Within this approach, one could use Raman scat- 
tering, which has a much better energy and momentum 
resolution than x-ray and neutron scattering. This ap- 
proach is feasible for graphene, which has a KA for the 
Raman-active E2g F-phonon [5] (Raman G-band). 

In this paper, we compute the variation of phonon 
frequency of the Raman G-band (E2g mode at T) 



in a graphene monolayer, as a function of the Fermi 
level. First, the calculations are done using a fully ab- 
initio approach within the customary adiabatic Born- 
Oppcnheimer approximation. Then, time-dependent 
perturbation theory (TDPT) is used to go beyond. 

Ab-initio calculations are done within density func- 
tional theory (DFT), using the functional of Ref. 
plane waves (30 Ry cutoff) and pseudopotentials 
The Brillouin zone (BZ) integration is done on a uni- 
form 64 X 64 X 1 grid. An electronic smearing of 0.01 Ry 
with the Fermi-Dirac distribution is used [8|. The two- 
dimensional graphene crystal is simulated using a super- 
cell geometry with an interlayer spacing of 7.5 A (if 
not otherwise stated). Phonon frequencies are calcu- 
lated within the approach of Ref. [§], using the PWSCF 
code [loj . The Fermi-energy shift is simulated by consid- 
ering an excess electronic charge which is compensated 
by a uniformly charged back-ground. 

The dependence of the Fermi energy ep on the surface 
electron- concentration a is determined by DFT (Fig[T]). 
In graphene, the gap is zero only for the two equivalent 
K and K' BZ-points and the electron energy e can be 
approximated as e^. /7r(K -I- k) = ±/3fc for the tt* and tt 
bands, where k is a small vector. Within this approxi- 
mation, at T = K temperature 



sign(eF)^ 



sign(ei.)— £^ 10.36 lO^^cm^^ n\ 
eV^ 



where (3 — 5.52 eV-A from DFT, sign(x) is the sign of 
X and ep — at the tt bands crossing. We remark that, 
from Fig [1] the typical electron-concentration obtained 
in experiments [l|,l4] corresponds to an important Fermi- 
level shift (~0.5 eV). For such shift, the linearized bands 
are still a good approximation (Fig. [T]). 

The dependence of the graphene lattice-spacing a on 
cr, a(cr), is obtained by minimizing -F(cr, A) = [E{a, A) — 
E{0, Ao)]/A with respect to A, where E{a,A) is the en- 
ergy of the graphene unit-cell, A is unit-cell area and 
Aq = 5.29 A^ is the equihbrium A [Tli at zero a. E{a, A) 
is computed by DFT letting the inter-layer spacing, L, 
tend to infinity in order to eliminate the spurious interac- 
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FIG. 1: Graphene monolayer. Upper panel: ef as a function 
of the surface electron-concentration a from DFT calculations 
and from linearized bands (at T = 300 K). Lower panel: in- 
plane lattice spacing a as a function of a. The fitting function 
is Eq. [2] and the dashed line is from Ref. • 



tion between the background and the charged sheet Il2l . 
Aa((T) = [a(cT) — a(0)]/a(0) was determined in Ref. [13| 
for intecalated graphite on the basis of a semi empiri- 
cal model. Using the same functional dependence as in 
Ref. [3], our DFT calculations are fitted by 



Aa(cr) = 6.748 ■ 10 



-6|^|3/2 



-I- 1.64- 10-V, 



(2) 



With cr = 3 10^3 
the lattice spacing variation is ~ 0.05%, which is, as ex- 
pected, of the same order of the valence-charge variation. 
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FIG. 2: Frequency of the E2g F phonon (Raman G band) 
as a function of cr: shift with respect to the zero-doping fre- 
quency. Calculations are done using standard DFT (adia- 
batic) or TDPT (dynamic), keeping the lattice-spacing con- 
stant (constant lattice) or varying it according to Eq. [2] (ex- 
panded lattice). Points are DFT calculations. Dashed line 
is from Eq. |3l Experiments should be compared with the 
continuous line. 



The frequency of the Y^ig V phonon is computed by 
static perturbation theory of the DFT energy i.e. 
from the linearized forces acting on the atoms due to the 
static displacement of the other atoms from their equilib- 
rium positions. This approach is based on the adiabatic 
Born-Oppenheimer approximation, which is the standard 
textbook approach for phonon calculations and is always 
used, to our knowledge, in the ab-initio frequency calcu- 
lations. The computed zero-doping phonon frequency is 
Lo^/{2'iTc) = 1554 cm~^, where c is the speed of light. The 
frequency variation Aw with a is reported in Fig. [21 Cal- 
culations are done keeping the lattice-spacing constant 
at a(0), or varying it according to Eq. O In this latter 
case, Acij is fitted by 

— = -2.13cr-0.0360cr2-0.00329c^3-0.226|crp/^ (3) 
27rc 

where a is in 10^'^ cm^^ and Aci;/(27rc) is in cni^^ units. 
The lattice-parameter variation is important, since it 
nearly doubles the frequency shift. However, Fig. [2] 
does not show the sudden increase of the phonon fre- 
quency with \a\, expected from the displacement of the 
KA wavevector with the doping. In particular, for a = 
3 10^^ cm^2, the frequency variation is ~ —0.5%, which 
excludes a magnification effect related to the KA. 

It is important to understand whether the absence of 
the KA is an artifact of the adiabatic approximation, 
used so far. Thus, we consider that a phonon is not a 
static perturbation but a dynamic one, oscillating at the 
frequency uj, which can be treated within time- dependent 
perturbation theory. Using such dynamic approach in the 
context of DFT [14| , the dynamical matrix of a phonon 
with momentum q, projected on the phonon normal- 
coordinate is 

^^{lo) = F^^{uj)+ I n(r)AV(r)dV 



- j An;(r)X(r,r')Anq(r')dVdV, (4) 

where n(r) is the charge density, A^y^ is the second 
derivative of the bare (purely-ionic) potential with re- 
spect to the phonon displacement. An is the derivative 
of 71, if(r,r') = 5''EH.c[n]/{5n{v)5n{v')), Eh a;c[n] is the 



Hartree and exchange-correlation functional, and 

I [/(k+q)m /kn] 



\D 



(k+q)m,kn 



e(k+q)m - ekn + flh) + iS 



(5) 



Here a factor 2 accounts for spin degeneracy, the 
sum is performed on iVk wavevectors, £'(k+q)m,k)i = 
((k -I- q)m| Ay|kn) is the electron-phonon coupling 
(EPC), Ay is the derivative of the Kohn-Sham poten- 
tial, |k7i) is a Bloch eigenstate with wavevector k, band 
index n and energy ek„, /kn = /T(ekn - ef), where /t is 
the Fermi-Dirac distribution and (5 is a small real number. 
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Imposing w = and S = in Eq. |4l one obtains the 
standard adiabatic approximation [9| and the phonon fre- 
quency is — ■y/Pq^(0)/Af, where M is the atomic 
mass. In the dynamic case, lo has to be determined self- 
consistently from u! — ^ Vif iuj) / M . However, consid- 
ering dynamic and doping effects as perturbations, at 
the lowest order one can insert the adiabatic zero-doping 
phonon frequency wJJ in Eq[4]and obtain the real part of 
the dynamic frequency from uj'^/ — 5Re [w^ ) /M . 

Let us consider the q ^ limit in Eq. [S] In the adia- 
batic case 



In the adiabatic case 



2 |-Dkm,knP[/km ~ /kn] 

^k , fkm - Ckn 



kn,kn I <5T(ek« - e_F), (6) 



where 5t{x) — —dfrix)/ {dx). In the dynamic case 



JVk 



E 



l-^krn^kn | [/km /kn] 
Ckm - Ekn + tlhJ^a + ' 



(7) 



In Eq. [S] (adiabatic case), there are two contributions, 
the first from inter-band and the second from intra-band 
transitions (depending on 5t and proportional to the den- 
sity of states at ep). On the contrary, in Eq. [7] (dynamic 
case) only inter-band transitions contribute. 
The variation of oj'^'^ with ei? is 



Aw 



(8) 



where is assumed that Aw <C The presence of a 

Kohn anomaly is associated to a singularity in the elec- 
tron screening, which, within the present formalism, can 
occur if the denominator of Eq. [5] approaches zero, i.e. for 
electrons near the Fermi level. Let us call F'^^ {lo) the part 
of i^Q ^ (w) obtained by restricting the k-sum on a circle 
of radius k centered on K, with {(3k —\ep\~ hw^) :s> k^T. 
The anomalous Aw is obtained by substituting T) with 
# in Eq.[5] 



AWa 

Awd = 5Re 



F"^(0) -F°(0) 



2MwO 

>-^(wg)-F°(wg) 
2A/w" 



(9) 
(10) 



in the adiabatic (Awa) and dynamic (Awd) cases. 
An analytic expression for F is obtained by i) lin- 
earizing the band dispersion; ii) writing the EPC as 
l^(K+k)«,(K+k)mP = {Dl)[l ± cos(26')], where 9 is the 
angle between the phonon-polarization and k, the sign 
± depends on the transition (see Eq. 6 and note 24 of 
Ref. jil) and {Dl) = 45.6 {eVf/k-'^ from DFT [Bl; iii) 
substituting 1/A^kEk ^i^^ 2Ao/(27r)2/ d^k in Eqs.[MZ], 
a factor 2 counts K and K', and k is measured from K. 



fe 

F'^-(O) = a / kdk 



MPk - €f) - M-Pk - ep) 
k 

5T{l3k-ep)-5T{-pk-ep)}, (11) 



where a = 2Ao(-Dr)/7r. Substituting Eq. [TT] into Eq. [HI 
one obtains Awa. At any T, Awq = 0. This re- 
sult is not trivial and comes from the exact cancella- 
tion of the inter-band (tt to tt*, first line of Eq. [TT|) 
and intra-band (tt to tt and tt* to tt*, second line of 
Eq. [TT|) . For example, at T = 0, both contributions 
to hl^uja are large and equal to a'|ei?| and —a'\ep\^ re- 
spectively, where a' = ha/ {2Muj^^l3'^) = 4.43 10"^ and 
a' /{2'iThc) — 35.8 cm~^/(ey). Concluding, an adiabatic 
calculation of w'^^ does not show any singular behavior 
in €p related to the Kohn anomaly, in agreement with 
the statc-of-thc-art adiabatic DFT calculations of Fig. [51 
In the dynamic case 

F^^(u^\ ^at /T(/3fc-^F)-/T(-/?fc-6^) | I 
Substituting Eq. [Hinto Eq. [TUl for T = 0, 



(12) 



-^In 



2 



~2~ 



(13) 



In this case, the situation is very different since the large 
inter-band contribution is not canceled by an intra-band 
term. In particular, there are two logarithmic divergences 
for ep = ±/iw"/2 and for \ep\ 3> ^a/2 the frequency 
increases as Q;'|ei?|. 

Awrf computed in this way takes into account tran- 
sitions between states close to the Fermi level. How- 
ever, the frequency is also affected by the variation of 
the lattice-spacing, by the transitions involving a state 
far from e^;' and by the second and third terms in Eq. U] 
All these contributions are accurately described by our 
adiabatic DFT calculations. Therefore, to compare with 
experiments, we add Aw^ to the adiabatic DFT fre- 
quency shift of Eq. [31 The results are shown in Fig. [51 for 
T = 300 K, and in Fig.[3|as a function of T for a smaller 
a range. Even at room temperature, the non-adiabatic 
Kohn anomaly magnifies the effect of the doping and for a 
valence-charge variation of -0.2% (-1-0.2%), the frequency 
varies by -1-1.5% (-1-0.7%). Aw is asymmetric with respect 
to ep and has a maximum for a ~ +3.5 10^^ cm~^. Since 
Awd is an even function of ep, this lack of electron-hole 
symmetry is entirely due to the adiabatic DFT contribu- 
tion. The ep = ±Sw°/2 logarithmic anomalies are visible 
at r = 4 and 70 K. The presence of a logarithmic KA 
in this two-dimensional system is quite remarkable since 
such divergences are typical of one-dimensional systems. 
They are present in graphene because of its particular 
massless Dirac-like electron band dispersion. 
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Finally, the Raman G-band has a finite homogeneous 
linewidth due to the decay of the phonon into electron- 
hole pairs. Such EPC broadening can be obtained either 
from the imaginary part of the TDPT dynamical matrix 
fEg. fT^ or, equivalently, from the Fermi golden rule [l^: 

(14) 

where 7 is the full-width half-maximum (FWHM) in 
cm~^. At T = and ep = 0, one recovers the result of 
Ref. [iBl, 7 — 11-0 cm^^. The phonon-phonon scattering 
contribution to the FWHM is smaller ('^1 cm~-^ [l3|) and 
independent oi ep- The total homogeneous FWHM is re- 
ported in Fig. [3l The FWHM displays a strong doping 
dependence; it suddenly drops for \a\ ~ 0.1 10^'^ cm^^ 
(|e_F| ~ 0.1 eV). Indeed, because of the energy and mo- 
mentum conservation, a T phonon decays into one elec- 
tron (hole) with energy fta;"/2 above (below) the level 
crossing. At T = QK such process is compatible with the 
Pauli exclusion-principle only if IefI < 
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FIG. 3: (Color online) Linewidth and dynamic frequency of 
the E29 r mode (Raman G band). See the caption of Fig. [1] 

Concluding, a Kohn anomaly dictates the dependence 
of the highest optical-phonon on the wavevector q, in 
undoped graphene Q. Here, we studied the impact of 
such anomaly on the q = phonon, as a function of the 
charge-doping a. We computed, from first-principles, the 
phonon frequency and linewidth of the E2g, F phonon 
(Raman G band) in the cr-range reached by recent FET 
experiments. Calculations are done using i) the custom- 
ary adiabatic Born-Oppenheimer approximation and ii) 
time-dependent perturbation theory to explore dynamic 
effects beyond this approximation. The two approaches 
provide very different results. The adiabatic phonon fre- 
quency displays a smooth dependence on a and it is not 
affected by the Kohn anomaly. On the contrary, when dy- 
namic effects are included, the phonon frequency and life- 



time display a strong dependence on a, due to the Kohn 
anomaly. The variation of the Raman G-band with the 
doping in a graphene-FET has been recently measured 
by two groups [isl . [l^ . Both experiments are well de- 
scribed by our dynamic calculation but not by the more 
approximate adiabatic one. We remark that the adiabatic 
Born-Oppenheimer approximation is considered valid in 
most materials and is commonly used for phonon calcu- 
lations. Here, we have shown that doped graphene is a 
spectacular example where this approximation miserably 
fails. 

We aknowledge useful discussions with A.M. Saitta, 
A.C. Ferrari and S. Piscanec. Calculation were done at 
IDRIS (Orsay, France), project n° 061202. 
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